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Abstract 

I review recent work on the "phase diagram" of the one-dimensional 
complex Ginzburg-Landau equation for system sizes at which chaos is 
extensive. Particular attention is paid to a detailed description of the 
spatiotemporally disordered regimes encountered. The nature of the 
transition lines separating these phases is discussed, and preliminary 
results are presented which aim at evaluating the phase diagram in 
the infinite-size, infinite-time, thermodynamic limit. 
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1. A prototype for studying spatiotemporal chaos 

The complex Ginzburg-Landau equation (CGL): 

d t A = A + (1 + ih)AA - (6 3 - i)\A\ 2 A 

where A is a complex field, x G [0, L] and 63 > 0, is of considerable importance 
to everybody interested in spatially extended non-equilibrium systems 5 . It 
accounts for the slow modulations, in space and time, of the oscillatory state 
in a physical system which has undergone a Hopf bifurcation 10 . As such, the 
CGL is closely related to numerous experimental situations 5 . The resulting 
universality and genericity are accompanied by specific features which make 
this model interesting for its own sake, simply as a prototype of spatially 
extended dynamical systems. 

There exist two important limits to CGL: the dissipative, relaxational 
"real Ginzburg-Landau" equation {b\ = 0, 63 — > 00) and the dispersive, in- 
tegrable nonlinear Schrodinger equation (NLS) (63 = 0, &i — > 00). Studying 
CGL in a comprehensive manner, one encounters dynamical regimes where 
the relative importance of dissipation and dispersion can be tuned at will in 
a non-integrable, non- variational system. This has been done recently for the 
one-dimensional case, uncovering the "phase diagram" presented in figure 1 
and discussed below 12 ' 2 . 

The relative simplicity of figure 1 stems mostly from the existence of a 
"thermodynamic limit" for the spatiotemporally chaotic regimes observed. 
As a matter of fact, away from the intricacy of the bifurcation diagrams 
at small sizes (L^,50), there exists a large-size limit beyond which chaos 
becomes extensive and can be characterized by intensive quantities indepen- 
dent of system size, boundary conditions, and, to a large extent, initial con- 
ditions. Perhaps the most convincing illustration, for CGL, of this essential 
property of spatiotemporal chaos can be found in a recent paper by Egolf and 
Greenside 6 where they show that the Lyapunov dimension is proportional to 
the system size L. 

Here, I review the various disordered "phases" observed in the one-dimen- 
sional case in the large-size limit and the nature of the transitions leading to 
them. Emphasis is put on the physical-space structures and objects which 
compose the spatiotemporal chaos as they could well play a crucial role when 
trying to build a statistical analysis of the disordered phases. This should 
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also hopefully provide a clearer picture of the "elementary mechanics" at 
play in these regimes. 

All the results reported and presented here were obtained using a pseudo- 
spectral code with periodic boundary conditions and second-order accuracy 
in space and time. Spatial resolution was typically 1024 modes for a domain 
of size L = 1000. Time steps varied like 1/y/bH, with typical values as 
indicated in the figure captions. 

2. Phases 

Early work on CGL has dealt with the problem of the linear stability of its 
family of plane-wave solutions A = a k expi(kx + Ukt) with a\ — (1 — k 2 )/b 3 
and tot = I/63 — (61 + 1 /6 3 )fc 2 . All these solutions are unstable for b\ > 63, a 
condition which defines the so-called "Benjamin- Feir line" (BF). For b\ < b 3 , 
plane-wave solutions with k 2 < &| ckhaus = (63 — b 1 )/(3b 3 — &i + 2/6 3 ) are 
linearly stable 8 . The BF line thus provides a first separation of the (61,63) 
plane, indicating where the constant-modulus, phase- winding solutions might 
play a role. 

I now give a brief account of what is known about the disordered phases. 
For all these spatiotemporal chaos regimes, chaos is extensive: the Lyapunov 
dimension (and the other quantities derived from the Lyapunov spectrum 
such as the entropy or the number of positive exponents) scales linearly with 
the system size L. 6 

2.1 Defect turbulence 

Above the BF line and to the left of Li, a strongly disordered phase is ob- 
served, best defined by a finite density of space-time zeros of \A\. This 
"defect turbulence", also named "amplitude turbulence", is characterized 
by a quasi-exponential decay of space and time correlation functions, and 
slower-than-exponential tails of the pdf of the phase gradient 0a;- 12 

One elementary process is the key ingredient of this strong spatiotemporal 
chaos: pulses of \A\ grow under the effect of the dispersion term; this "self- 
focusing" is stopped by the action of dissipation, breaking the pulse. Such 
events can be seen in the spatiotemporal plots presented in figure 2; regions 



3 



of almost zero amplitude are present, another indication that pulses are the 
relevant objects to consider when approaching the NLS limit. As one gets 
closer to the NLS (increasing bi), the pulses get sharper and higher, whereas 
closer to the BF line, dissipation dominates and \A\ rarely overpasses the 
"saturation" value 1/^/%. 

This disordered phase is reached from (almost) all initial conditions for 
parameter values to the left of L 3 . In the triangular region delimited by BF, 
!_! and L 3 . it coexists with "phase turbulence", which is described below. 

2.2 Phase turbulence 

Phase turbulence is a weakly disordered regime observed in the region of 
parameter space above the BF line to the right of l_i, as well as in the 
triangular region delimited by BF, l_i and L3. It is best defined by the 
absence of space-time defects or, equivalently, by the bounded character of 
the pdf of 4> x and \A\ (see below for a discussion of this point). The absence 
of phase singularities implies that the "winding number" v = / L <f) x dx is a 
conserved quantity equal to a multiple of 2n which classifies the different 
attractors reached at given values of b\ and b 3 . 

Chaos is very weak, as shown by the slow decay of space and time cor- 
relation functions, indicative of diffusive or sub-diffusive modes. 12 This is 
corroborated by the flat shoulder of the spatial power spectrum of \A\ at 
low wavenumbers, reminiscent of the Kuramoto-Sivashinsky equation (KS). 
The dynamics is in fact very similar to that of KS, which is not surprising 
since this equation was originally derived to describe the phase dynamics of 
CGL near the BF line. 10 The close relationship with KS suggests, in turn, 
that the Kardar-Parisi-Zhang interface equation (KPZ) 9 , which have been 
argued to account for the large-scale, long-time, properties of KS, 13 might 
provide the correct description of the long- wavelength limit of phase turbu- 
lence. In such a setting, one expects exponential decay of spatial correlations 
and stretched-exponential decay of temporal correlations. 9 

Spatiotemporal diagrams (figure 3) reveal that the objects involved are 
"shocks" of \A\ (similar graphs for <f> x show localized modulations of the 
wavenumber, and look essentially the same). Close to the BF line, the dy- 
namics is indeed very much like in KS (figure 3a), while away from it trains 
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of propagative shocks become more and more frequent (figure 3b). The over- 
all effect of v is to introduce a general drift on the localized shocks of \A\ 
which increases with abs(z/) (figure 3c). The phase turbulence regimes thus 
look very different depending on v and the distance, in the parameter plane, 
from (bi, 6 3 ) to the BF line. Increasing this distance, the characteristic scales 
decrease (the difference of the time-scales between figures 3a and 3b is only 
partly due to the variation of the basic frequency like y^). The other impor- 
tant effect, the appearance of trains of propagative shocks, is an indication 
that KS is not a valid description away from the BF line and that other 
terms, incorporating odd derivatives of 0, are present in the "effective phase 
equation" describing the corresponding phase turbulence regimes. Whether 
the connection to KPZ still holds in such circumstances and what is then the 
valid phase equation replacing KS, are open problems under investigation. 

It must be mentioned, at this point, that it is still a controversial and 
unresolved issue to know whether phase turbulence "exists" in the thermo- 
dynamic limit, i.e., whether it is a transient phenomenon or a finite-size 
artifact. In sections 3.1 and 3.3 below, I discuss this point in relationship 
with the nature of the transition lines delimiting the phase turbulence region 
in figure 1. In any case, even if there always exists a very small (essentially 
undetectable) density of defects, the spacetime regions between those rare 
events are large enough so as to justify the study of phase turbulence for its 
own sake. 

2.3 Spatiotemporal intermittency 

Below the BF line, plane- wave solutions of wavenumber k 2 < &l ckhaus are 
linearly stable. This does not preclude the existence, to the left of L 2 , of 
other, mostly chaotic, solutions which are easily reached for initial conditions 
outside the basin of attraction of the stable plane waves. 2 These disordered 
regimes are spatiotemporal intermittency regimes: they consist of space-time 
regions of stable plane waves separated by localized objects evolving and in- 
teracting in a complex manner (figure 4). The plane waves constitute the 
passive, "absorbing" state while the localized objects carry the spatiotem- 
poral disorder. The simplest of these objects are members of the family of 
exact solutions found by Nozaki and Bekki. 1 These propagating "holes" of 
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\A\ are not defects, strictly speaking, even though \A\ might remain very- 
small in their core. 2 Defects (zeros of A) do occur, but at scattered points 
in spacetime. This makes the density of defects difficult to estimate and 
rather worthless, since, for example, the average number (in space) of hole- 
like objects is probably a better characterisation of the dynamics. Likewise, 
space and time correlation functions are not the easiest quantities to measure 
when evaluating the coherence scales in these regimes, due to the intermit- 
tent character of the disorder. As usual with spatiotemporal intermittency, 4 
one can take advantage of the intrinsic binary structure of the spatiotempo- 
ral dynamics and evaluate, say, the distributions of sizes and lifetimes of the 
patches of plane wave solutions. These distributions are roughly exponential, 
yielding characteristic coherence scales. 2 

The spatiotemporal intermittency regimes do exhibit space-time defects. 
To that extent, they may be considered as defect turbulence regimes. The key 
difference with the region of defect turbulence described above is that here 
the defects do not appear spontaneously, they are produced by the localized 
objects carrying the disorder. 

2.4 The "bichaos" region 

The domain of parameter space delimited by the BF, l_i and L 3 lines deserves 
further comment. Depending on the initial conditions, one can reach one of 
two spatiotemporally chaotic regimes, phase turbulence or defect turbulence. 
A closer look at the defect turbulence reveals that it consists of localized 
propagating and branching hole-like objects which separate more quiescent 
space-time regions (figure 5). These objects do not appear spontaneously. 
The quiescent regions are nothing else than patches of phase turbulence, with 
their characteristic shocks of \A\. Defect turbulence is thus a spatiotemporal 
intermittency regime, with the absorbing state being the phase turbulence 
regime. 

3. Transitions 

I now examine the various transition lines between the phases and discuss to 
what extent they can be considered as phase transitions. 
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3.1 The Li line 



Line l_i separates defect turbulence from phase turbulence for values of 
bi ^1.85. As such, its best definition is given by the parameter values for 
which the density of defects no goes to zero. It was shown 12 that, for bi = 3.5, 
no ~ (b 3 — 6g) 2 with 63 ~ 1.29. No hysteresis has been detected. This tran- 
sition thus appears continuous, even though one cannot exclude a crossover 
scenario where would remain small beyond 63. 

The situation is not as clear from the point of view of the correlation 
length £ deducted from the exponential decay of the two-point correlation 
function: approaching l_i, £ increases but does not diverge at the transition 
point. 12 ' 7 Although difficult to measure in the phase turbulence regimes, in- 
direct arguments imply that £ remains finite. Only a change of behavior is 
observed, apparently accompanied by a discontinuity of d£/dfe 3 . A similar 
situation holds for the variation of the Lyapunov dimension and the other 
quantities derived from the Lyapunov spectrum when crossing l_i, as shown 
by Egolf and Greenside. 6 In other words, no critical behavior seems to be 
present at the transition, even though it is continuous. 

This remark is valid for the experimental conditions within which the 
above results were obtained. In particular, the system sizes used ensure the 
extensivity of chaos, but no assessment of finite-size effects has been made 
in order to evaluate the infinite-size, infinite-time, "thermodynamic" limit. 
The main question in this respect is that of the position of Lx as L — > 00. If 
l_i remains away from the BF line, the existence of phase turbulence in the 
thermodynamic limit is assured and the transition at l_i does not show any 
critical behavior. Another hypothesis consists in l_i moving toward BF as 
L increases, with the correlation length at threshold increasing and finally 
diverging for an infinite-size system. In this case, the phase turbulence region 
disappears in the thermodynamic limit and the BF line is the asymptotic 
boundary of defect turbulence. 

Preliminary results (figure 6) indicate a clear displacement of l_i toward 
BF as L increases, making the second hypothesis more likely. A similar 
trend was recently found by Egolf and Greenside. 7 Extensive numerical work 
is under completion to obtain a quantitative extrapolation of the position of 
l_i in the infinite-size limit. 
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3.2 The L 2 line 



The line l_2, best defined as the limit of existence of spatiotemporal disor- 
dered states, is difficult to determine in a precise manner. 2 This is mostly due 
to the appearance, in the transition region, of "frozen" states, i.e. spatially- 
disordered arrangements of localized objects with trivial time dependence 
(figure 7). The coherence scales of the spatiotemporal intermittency regimes 
increase when approaching l_2, but their divergence is not observed; the sys- 
tem "falls" on a frozen state. As often in spatiotemporal intermittency 
situations 3 , the deterministic features of the system "mask" the directed- 
percolation-like phase transition. The line l_2 is thus determined only crudely, 
awaiting further progress. In particular, two directions seem worth investi- 
gating: a "local" approach based on the study of the interactions between the 
localized objects at the origin of the spatiotemporal disorder, and a "global" 
approach based on Lyapunov analysis of the disordered regimes. 

3.3 The L 3 line 

The L 3 line delimits the parameter space region where phase turbulence 
ceases to coexist with defect turbulence. Here I give a preliminary account 
of ongoing work on the nature of the breakdown of phase turbulence occurring 
when crossing L 3 (say, by decreasing 63), for a winding number v — 0. 

The transition is hysteretic: defect turbulence persists when crossing 
L 3 back, up to line l_2. The breakdown of phase turbulence 11 occurs via 
the nucleation of a first defect, followed by the quasi-deterministic invasion 
of the phase turbulent state by the defect phase (figure 8). The speed of in- 
vasion is always finite and rather large; it is given by the average propagation 
speed vd of the localized objects composing the defect phase in this region of 
parameter space. The precursor of the first defect is apparently a local event: 
one of the shocks of \A\ composing the phase turbulence accelerates, and a 
depression of \A\ develops at its tail. When the velocity reaches v D) the first 
hole-like object typical of the defect phase is effectively created, and quickly 
generates the first space-time defect. It is difficult to determine what triggers 
the initial acceleration of one of the shocks of \ A\. As mentioned above, more 
and more trains of propagating shocks are observed in the phase turbulence 
regimes as one goes away, in parameter space, from the BF line. Simul- 



8 



taneously, the average velocity vs of these propagative shocks increases. A 
possible interpretation of the nucleation is that, for parameter values beyond 
L 3 , the fluctuations around vs bring the speed of one shock past a critical 
value, after which the shock is attracted to the hole-like solution characteris- 
tic of the defect phase. Such a critical value would be related to correponding 
critical values of the local minimum of \A\ and local extremum of <f> x , since as 
the velocity of the shock increases, the local depression of \A\ at its tail deep- 
ens, and the local phase gradient increases. Ongoing work aims at making 
those mostly qualitative statements more quantitative. 

As described above, the breakdown of phase turbulence strongly resem- 
bles a first-order phase transition. So far, L 3 has been determined numeri- 
cally by merely checking whether the breakdown occurs within a fixed (very 
long) integration time. 12 However crude, this methodology produced a sharp 
transition line, all the more so as no size effects could be detected. 

A quantitative estimate of the probability of nucleation is needed to im- 
prove this situation. Preliminary results on the study of the pdfs of \A\ 
and 4> x in the phase turbulence regimes provide the beginning of an answer 
(figure 9). Away from L 3 (figures 9a and 9c), these distributions have ei- 
ther Gaussian or strictly bounded tails. Close to L 3 (figures 9b and 9d), 
these tails are exponential, signalling a clear change of behavior on which 
the definition of L 3 can be firmly based. Nevertheless, the question of the 
Gaussian/bounded character of the tails away from L 3 remains a crucial 
one since it determines whether the probability of nucleation can be strictly 
zero or not. A plausible scenario would consist of exponential tails truncated 
by a finite bound moving to infinity (resp. zero) for <f> x (resp. \A\) when 
approaching L 3 . Whereas this scenario would assure the existence of the 
phase turbulence in the infinite-size, infinite-time limit, the other scenario of 
a crossover from Gaussian to exponential tails would not. Current numerical 
efforts aim at providing data to decide this. 

In the same spirit, another important question deals with the estimation 
of size-effects in this region of parameter space. So far, results for sizes L = 
1000 and L = 4000 have not shown any significant difference, strengthening 
the case for a first-order transition. 

Finally, the overall effect of a non-zero winding number v is to give rise 
to earlier breakdown of the phase turbulence regime, moving L 3 to the right, 
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as expected from the general additional drift of the shocks in this case. 
4. Open problems 

The "phase diagram" presented in figure 1 was obtained for systems sizes 
L large (for which the extensivity of chaos is ensured) but finite. From the 
above discussion of lines L x and L 3 , it appears that an estimation of the 
finite-size effects is necessary, and even crucial, since the existence of phase 
turbulence depends on the position and nature of these transition lines in 
the infinite-size, infinite-time, "thermodynamic" limit. For line L3, a de- 
tailed study of the pdfs of interest in the phase turbulence regimes should be 
able to provide an objective criterion which can help to determine precisely 
the breakdown of phase turbulence. In this respect, preliminary results (fig- 
ure 9) are encouraging. For the line l_2 delimiting the domain of existence of 
spatiotemporal intermittency, the key point lies in finding a quantifier which 
can be used to determine its position efficiently. Quantities related to the 
Lyapunov spectrum might constitute such a quantifier. 

Apart from these "global", statistical quantities, approaches based on 
the study of the local objects involved in the spatiotemporally disordered 
dynamics ("shocks" for phase turbulence, "pulses" and "defects" for defect 
turbulence, "sources" and "sinks" for spatiotemporal intermittency) might 
also help to clarify some of the problems left open. 

At any rate, and whatever the situation is in the thermodynamic limit, the 
"finite-size" phase diagram of figure 1 is of interest to all "experimentalists" 
(i.e., anybody working with finite-size systems observable during a finite 
time). In this respect, it shows how varied and complex statiotemporally 
chaotic regimes can emerge from such a simple equation as CGL. I hope the 
above review will also provide guidelines for dealing with similar cases, in 
particular when the question of what to measure arises. 

I thank D.A. Egolf and H.S. Greenside for stimulating exchanges and 
the communication of their unpublished results, and P. Manneville for many 
discussions and comments. 



10 



References 



1 Bekki, N., and K. Nozaki, "Exact Solutions of the Generalized Ginzburg- 
Landau Equation." J. Phys. Soc. Jap. 53 (1984): 1581-1582; "Forma- 
tion of Spatial Patterns and Holes in the Generalized Ginzburg-Landau 
Equation." Phys. Lett. A 110 (1985): 133-135. 

2 Chate, H., "Spatiotemporal Intermittency Regimes of the One-Dimensional 
Complex Ginzburg-Landau Equation", to appear in Nonlinearity 1993. 

3 Chate, H., and P. Manneville, "Role of Defects in the Transition to 
Turbulence via Spatiotemporal Intermittency." Physica D 37 (1989): 
33-41. 

4 Chate, H., and P. Manneville "Spatiotemporal Intermittency in Cou- 
pled Map Lattices." Physica D 32 (1988): 409-422; Chate, H., "Sub- 
critical Bifurcations and Spatiotemporal Intermittency." In: Sponta- 
neous Formation of Space- Time Structures and Criticality, edited by 
T. Riste and D. Sherrington, 273-311. Dordrecht: Kluwer, 1991. 

5 Cross, M.C., and Hohenberg, P.C., "Pattern Formation Outside of 
Equilibrium", to appear in Rev. Mod. Phys. 1993, and references 
therein. 

6 Egolf, D.A. and H.S. Greenside, "Complexity Versus Disorder for Spa- 
tiotemporal Chaos." Preprint available from the LANL database "chao- 
dyn@xyz.lanl.gov" under number # 9307010. See also the abstract of 
a poster by the same authors in this volume. 

7 Greenside, H.S., private communication. 

8 Janiaud, B., A. Pumir, D. Bensimon, V. Croquette, H. Richter, and 
L. Kramer, "The Eckhaus Instability for Travelling Waves." Physica D 
55 (1992): 259-269. 

9 Kardar, M., G. Parisi and Y.C. Zhang, "Dynamic Scaling of Growing 
Interfaces". Phys. Rev. Lett. 56 (1986) 889-892. 



11 



Kuramoto, Y., Chemical Oscillations, Waves and Turbulence. Tokyo: 
Springer, 1984; Lega, J., "Defauts Topologiques Associes a la Brisure 
de l'lnvariance de Translation dans le Temps." These de Doctorat, Uni- 
versite de Nice, 1989; Coullet, P., L. Gil, and J. Lega "Defect Mediated 
Turbulence." Phys. Rev. Lett. 62 (1989): 1619-1622. 

Sakaguchi, H., "Breakdown of the Phase Dynamics." Prog. Theor. 
Phys. 84 (1990): 792-800. 

Shraiman, B.I., A. Pumir, W. van Saarloos, P.C. Hohenberg, H. Chate, 
and M. Holen, "Spatiotemporal Chaos in the One-Dimensional Com- 
plex Ginzburg-Landau Equation." Physica D 57 (1992): 241-248; 
Pumir, A., B.I. Shraiman, W. van Saarloos, P.C. Hohenberg, H. Chate, 
and M. Holen, "Phase vs. Defect Turbulence in the One-Dimensional 
Complex Ginzburg-Landau Equation." In: Ordered and Turbulent Pat- 
terns in Taylor- Couette Flows, edited by CD. Andereck. New York: 
Plenum Press, 1992. 

Yakhot, V., "Large-Scale Properties of Unstable Systems Governed by 
the Kuramoto-Sivashinsky Equation." Phys. Rev. A. 24 (1981): 642- 
644; Zaleski, S., "A Stochastic Model for the Large Scale Dynamics of 
Some Fluctuating Interfaces." Physica D 34 (1989): 427-438; Procac- 
cia, L, M.H. Jensen, V.S. L'vov, K. Sneppen and R. Zeitak, "Surface 
Roughening and the Long- Wavelength Properties of the Kuramoto- 
Sivashinsky Equation." Phys. Rev. A. 46 (1992): 3220-3224; Sneppen, 
K., J. Krug, M.H. Jensen, C. Jayaprakash and T. Bohr, "Dynamic Scal- 
ing and Crossover Analysis for the Kuramoto-Sivashinsky Equation." 
Phys. Rev. A 46 (1992): 7351-7354. 



12 



Figure Captions 



Figure 1: Phase diagram in the (61,63) parameter plane. Lines U and L 3 were 
determined in 12 , line L 2 in 2 . See 2 for numerical details. 

Figure 2: Spatiotemporal representation of a defect turbulence regime in a system 
of size L = 250 for b\ = 3.5 and 63 = 0.5 (periodic boundary conditions, 
timestep St = 0.02). Time is running upward for AT = 82.5 (transient 
discarded). Part (a) shows the evolution of \A\ (grey scale between 
\A\ = 0/black and \A\ = 2/white). Part (b) shows the corresponding 
evolution of = arg(A) (grey scale between 0/white and 7r/black). 

Figure 3: Spatiotemporal evolution of \A\ in phase turbulence regimes in a system 
of size L = 1000 for bi = 3.5 (periodic boundary conditions, timestep 
St = 0.16, transient discarded, time running upward, grey scale with 
white corresponding to maximum value of \A\). Part (a): close to the 
BF line (y = 0, 63 = 2, AT = 3300, 0.68 < \A\ < 0.75). Part (b): 
close to Li {y = 0, b 3 = 1.4, AT = 1650, 0.69 < \A\ < 1.01). Part (c): 
same as (a) but with v = — 10n. 

Figure 4: Spatiotemporal intermittency regime in a system of size L = 250 (peri- 
odic boundary conditions, timestep St = 0.04) for 61 = and 63 = 0.5. 
Time is running upward for AT = 165 (transient discarded). Same 
representation as in figure 2, but < \A\ < 1.413 in part (a). 

Figure 5: Defect turbulence in the bichaos region (L = 1000, b\ = 1.5, 63 = 1, 
AT = 1312, periodic boundary conditions, timestep St = 0.04, and 
transient discarded). Spatiotemporal representation of | A\ in grey scale 
from \A\ = 1.12 (white) to \A\ = (black); time is running upward. 
Note that the hole-like propagating and branching objects typical of 
the defect phase do not appear spontaneously in the phase turbulent 
medium. 

Figure 6: Variation of the defect density n D (xlO 5 on the graph) with 6 3 near 
Li. These results were obtained from runs of duration AT ~ 20000 
after transients on systems of size L = 1000 and 2000 with periodic 
boundary conditions and timestep St = 0.04. Although these results are 
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still not very precise (at such low densities, much longer runs would be 
necessary), the general trend is significant: higher densities are reached 
for larger sizes. 

Figure 7: Frozen state observed near l_2 in a system of size L = 500 with b\ = 
—0.9 and 63 = 0.18 (periodic boundary conditions). This spatially 
disordered state is made of zero- velocity Nozaki-Bekki holes and shocks. 
\A\ is stationary, while winds regularly along time. 

Figure 8: Breakdown of phase turbulence past the L 3 line in a system of size L = 
1000 for bi = 1, b 3 = 0.675 and v — (periodic boundary conditions, 
timestep St = 0.16). Initial condition was a phase turbulent state for 
63 = 0.7. The breakdown occurred near t = 12000 as shown by the 
time series of \A\ (part (a)). Part (b) is the spatiotemporal evolution 
of \A\ for t G [9600, 12800] (grey scale from \A\ = 1.29/white to \A\ = 
1.1 /black, time running upward); at this resolution, the details of the 
defect phase are not visible. Part (c): zoom of (b) with a grey scale 
and a resolution adapted to the defect phase (x G [500,1000], t G 
[12550,12800], < \A\ < 1.29); the nucleation of the first defect is 
obvious. Part (d): as in (b) but for the locus of the minimum (in 
space) of \A\. 

Figure 9: Probability distribution functions of \A\ for phase turbulence regimes 
in the bichaos region. These distributions were obtained from the sim- 
ulation of a system of size L = 4000 during AT = 11500 sampled every 
4 timesteps St = 0.16. (periodic boundary conditions, transient dis- 
carded). Parts (a) and (c): v — 0, b± — 1, b 3 — 0.9. Parts (b) and (d): 
v = 0, 61 = 1, 63 = 0.7. 
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